A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

11.06.2; 11.09.1 Phoenix; 11.12.1; 11.19.5; 11.19.6 



ASTRONOMY 

AND 
ASTROPHYSICS 



Stellar populations in the Phoenix dwarf galaxy. 



Enrico V. Held^, Ivo Saviane^, and Yazan Momany^ 

^ Osservatorio Astronomico di Padova, vicolo deH'Osservatorio 5, 1-35122 Padova, Italy 

^ Universita di Padova, Dipartimento di Astronomia, vicolo deH'Osservatorio 5, 1-35122 Padova, Italy 



> ■ 

i> : 
o 

m 
o 

ON 
Oh 

6 ■ 

■ 



X 



Abstract. We have obtained deep CCD photometry in the B, 
V, and / bands of Phoenix, a galaxy considered a transition 
case between dwarf spheroidal (dSph) and dwarf irregular (dl) 
galaxies. A comparison of our data with the giant branches 
of Galactic globular clusters gives a mean metal abundance 
[Fe/H]= —1.81 ± 0.10 dex. The presence of an intrinsic color 
dispersion in the upper red giant branch (RGB) suggests an 
abundance range of about 0.5 dex, although a range in age 
may also affect the RGB width. The color-magnitude diagram 
(CMD) of Phoenix reveals for the first time a horizontal branch 
(HB) predominantly red yet moderately extended to the blue, 
similar to those of Leo II or And I, at V « 23.8. The detection 
of a relatively blue HB indicates the presence of a significant 
population with age comparable to that of old halo Galactic 
globular clusters. As in other dwarf spheroidals, this HB mor- 
phology in a metal-poor system indicates a mild "second pa- 
rameter" effect. The mean level of the HB has been used to 
derive a true distance modulus 23.21 ± 0.08, in good agree- 
ment with the distance modulus 23.04 ± 0.07 estimated from 
the well defined cutoff of the red giant branch at / « 23.1. 
This confirms the correct identification of the RGB tip. We 
also find a radial gradient in the Phoenix HB morphology, as 
measured by an increasing ratio of blue HB stars to red giant 
stars in the outskirts of the galaxy. The color-magnitude dia- 
grams show a small number of stars above the tip of the RGB, 
well in excess over field contamination, that most likely are 
asymptotic giant branch (AGB) stars belonging to an interme- 
diate age population. Their number indicates that the fraction 
of intermediate age (3 to 10 Gyr) population in Phoenix is ap- 
proximately 30^0%. A young stellar population is definitely 
present in Phoenix, consistent with a star formation episode 
started at least 0.6 Gyr ago, up to 1 x 10^ yr ago. Both young 
stars and AGB stars are centrally concentrated, which indicates 
that recent star formation preferentially occurred in the inner 
galaxy regions. In many respects, including an extended star 
formation history and even the presence of a modest amount of 
neutral hydrogen. Phoenix appears not dissimilar from dwarf 
spheroidal galaxies in the Local Group. 
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1. Introduction 

Photometric studies of dwarf spheroidals, and in particular of 
their color-magnitude diagrams, have shown a variety of star 
formation histories, ranging from the globular cluster-like dia- 



grams of Draco or Tucana (Grillmair et al. 1997; Saviane et al 



1996, hereafter Paper I) to the considerable intermediate age 



populations of Cari na (Sm ecker-Hane e t al. 1994) and Fornax 



(Beauchamp et al. |l995| ; Stetson et al. |l998| ; Saviane & Held 



19991 ), with a range of intermediate cases f or whi ch we refer 



t he rea der to the recent reviews of Da Costa ( |1998| ) and Mateo 
( |199S| ). 

The dwarf galaxy Phoenix represented a particularly inter- 
esting case in this respect since it is currently forming stars 
on top of a seemingly old, metal-poor population. Ortolani 
& Gratton ( |l988i hereafter OG88) obtained B, V CCD pho- 
tometry in a 2' X 3' field centered on the blue star associa- 



tion first described by Canterna & Flower (1977). Their color- 
magnitude diagrams show bright blue stars likely belonging to 
a 10^ yr-old recent burst of star formation, involving ~ 10^ 
A/0, superposed onto a predominantly old (> 10^° yr) metal- 
poor stellar population. A comparison with metal-poor Galac- 
tic globular clusters led OG88 to estimate a distance modulus 
(to — Al) w 23 — 23.5, thus locating Phoenix well within the 
Local Group. OG88 also found ?l{B — V) color distribution on 
the RGB somewhat broader than accounted for by instrumental 
eiTors, the intrinsic color dispersion corresponding to an abun- 
dance range of about 0.6 dex. A further study of the resolved 



stellar populations of Phoenix by van de Rydt et al. ( 1 99 1 here- 



after VDK91), yielded a slightly shorter distance modulus of 
(to — Af)o = 23.1 and a metal abundance [Fe/H]= —2.0. 
The integrated luminosity. My ~ —9.9, is similar to that of 
Carina. These photometric studies together with limited spec- 



troscopic data of Da Costa (1994), also indicated the presence 



of low-luminosity carbon stars consistent with an intermediate 
age population older than 8-10 Gyr. 
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Both studies concluded that Phoenix may belong to an in- 
termediate class between dl's and dSph's, a view supported by 
the detection of neutral hydrogen towards the optical galaxy 
at a heliocentric radial velocity of 56 km s^^ (Carignan et al. 



1991[ ). More recently. Young & Lo ( |1997[ ) obtained VLA maps 
of Hi in a region around Phoenix at a resolution ^ 2'. Their 
maps show an Hi cloud ("cloud B") at about 8' — 10' from the 
optical center, at a heliocentric velocity of 50 — 60 km s^^, 
whi ch is identified with the emission detected by Carignan et 
al. (1991). Its velocity dispersion is low, only ~ 3 km s^^. In 
addition, that study revealed a new Hi feature peaked at ^ —23 
km s^^, about 5' southwest of the optical galaxy ("cloud A"), 
whose curved shape suggests a physical link with the galaxy 
(see the Hi contour maps in their Fig. 17). However, it is con- 
ceivable that both clouds are associated with the complex Hi 
velocity structure of the Magellanic Stream. 

The presence of young stars along with neutral gas in an 
essentially old galaxy may provide the clue to understanding 
the role of multiple bursts of star formation in the evolution of 
dwarf spheroidals (and possibly of their more massive counter- 
part, the dwarf ellipticals). For this reason, we have undertaken 
a study of the old, young, and intermediate age populations of 
Phoenix. A first task was establishing whether Phoenix harbors 
a significant populations of metal-poor stars of age comparable 
to Galactic globular clusters, in which case a blue HB is ex- 
pected. Indeed, we have been able to detect for the first time a 
moderately blue HB atV « 23.7. A further aim was a quan- 
titative estimate of the contribution of young and intermediate 
age stars to the stellar populations of Phoenix, in particular to 
understand if the recent burst of star formation represents an 
isolated episode after a long quiescent period. 

The plan of the paper is as follows. Section ^ presents our 
new observations in B, V, and / of the Phoenix dwarf over 
a 6'.9 X 6'.9 field. Careful analysis of the instrumental errors 
and foreground/background contamination, essential to esti- 
mate the contribution of the young and intermediate age stellar 
populations, is also presented. In Sect. ^ we describe the main 
features of the color-magnitude diagrams, and outline the spa- 
tial distribution of stars of different ages. In particular, we in- 
troduce our first detection of the horizontal branch of Phoenix. 
The basic properties of the galaxy are derived in Sect. 0. There 
we provide a new distance estimate based both on the / magni- 
tude of the luminosity of the RGB tip, and on the mean V mag- 
nitude of horizontal branch stars. The wide photometric base- 
line employed in this study allows us to re-address the problem 
of an intrinsic color range and abundance dispersion of the red 
giants. The stellar populations of the Phoenix dwarf are dis- 
cussed in Sect. |^, where in particular we estimate the contribu- 
tion of young and intermediate age stars to the galaxy luminos- 
ity, and map the distribution of the star-forming regions across 
the galaxy. The main results of this paper are summarized in 
Sect. along with concluding remarks. 



2. Observations and data reduction 

Observations of a single 6. '9 x 6'.9 field centered on Phoenix 
were obtained on September 7 and 8, 1996 using EFOSC2 at 
the ESO/MPI 2.2m telescope. The detector was a thinned Lo- 
ral CCD with 2048 x 2048 pixels of 0'.'26. The images were 
trimmed to 1600 x 1600 pixels during readout to avoid the vi- 
gnetted edges of the image. Fast readout, appropriate for broad 
band imaging, yielded a read-out noise of 7.8 e^/pixel and a 
conversion factor of 1.34 e^/ADU. A control field located at a 
few degrees from the galaxy was also observed to estimate the 
foreground and galaxy background contamination. Sky condi- 
tions were photometric on the second night of the run. The see- 
ing was quite variable, so that only the best pairs of images in 
each filter were selected for processing. Table |l] lists the images 
used in this study, including the comparison field exposures. 
Images with poorer seeing (not listed here) were used for cali- 
bration purposes, namely to evaluate the zero point stability. 

Table 1. The journal of observations 



ID 


Night 


Filter 


icxp [s] 


airmass 


FWHM" 


Ph/b5 


2 


B 


2400 


1.04 


1.4 


Ph/b6 


2 


B 


2400 


1.04 


1.3 


Ph/i2 


2 


I 


2400 


1.06 


1.2 


Ph/i3 


2 


I 


1800 


1.10 


1.3 


Ph/vl 


1 


V 


1500 


1.29 


1.5 


Ph/v2 


1 


V 


1500 


1.21 


1.4 


bck/bl 


1 


B 


2400 


1.04 


1.3 


bck/b2 


1 


B 


2400 


1.06 


1.3 


bck/il 


1 


I 


2400 


1.07 


1.3 


bck/i2 


1 


I 


2400 


1.04 


1.3 


bck/vl 


1 


V 


1500 


1.15 


1.7 


bck/v2 


1 


V 


1500 


1.11 


1.7 



Image pre-processing was carried out with the ESO MIDAS 
package in the standard way. First, all images were cleaned us- 
ing a map of the bad features of the CCD. After checking the 
bias stability, a master bias was subtracted from all images. 
Several daylight dome flats and twilight flat-fields were ob- 
tained in each filter and used to construct dome and sky master 
flat-fields. The sky flats were preferred for correcting the sci- 
entific frames, because of the higher signal-to-noise ratio and 
better color match, though the dome flats were almost as good. 
The images were then registered to a common reference frame 
and coadded. The size of the Point Spread Function (PSF) of 
the averaged images shows no degradation with respect to the 
individual frames. As an example, our V sum image is shown 
in Fig. [l|. 

Stellar photometry was performed using DAOPHOT II and 



ALLSTAR (Stetson 1987). The PSF was iteratively constructed 
from a starting list of about 100 (relatively bright) reference 
stars, uniformly distributed across the CCD. After fitting a pre- 
liminary PSF to the sum images, all reference stars with faint 
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Fig. 1. A centi-al field (5.'3 x 5.'3) of the coadded V image of 
Phoenix. North is to the top and East to the left. The ou tlines 
indicate the regions studied by Ort olani & Gratton ( 1988 ; inner 
box) and van de Rydt et al. (1991 outer box) 



neighbors were culled out of the catalog by careful visual in- 
spection of the star-subtracted images. A better point spread 
function was then obtained from the selected list of PSF stars. 
Different PSF shapes were tested by checking the residuals 
after subtraction. The best fitting was obtained using a Mof- 
fat function with /3 = 1.5 and quadratic radial dependence. 
The subtraction residuals were generally satisfactory due to the 
good sampling of stellar images. Photometry was then obtained 
by running DAOPHOT and ALLSTAR twice on the sum images. 

The photometric errors and the degree of completeness of 
our data were estimated by extensive artificial star experiments. 
We used a technique in which artificial stars added to each im- 
age are randomly scattered about the vertices of a fixed grid or 
"lattice". This prevents overlapping of simulated stars ("self- 
crowding"). Lists of input stars (typically ~ 900) were cre- 
ated for the V image, with randomly distributed coordinates 
and V magnitudes uniformly distributed within the limits of 
the Phoenix CMD. Random {B — V) and (V — I) colors were 
employed, so that the artificial star colors spanned the entire 
useful intervals in the Phoenix color-magnitude diagrams. We 
made 30 experiments per filter, corresponding to ~ 27000 stars 
in each bandpass. The frames containing the simulated stars 
were reduced exactly in the same way as the real images. The 
master catalog of star experiments contains the artificial stars 
retrieved in at least one filter, ready for calibration as the real 
photometry catalogs. 

Table || lists the photometric errors obtained in each band- 
pass separately. The first column gives the center of the mag- 
nitude bin in B, V, and /. The photometric errors were ob- 



Table 2. The photometric errors (la) determined from the ar- 
tificial star experiments 



bin 


cry 






iO.Zj 








1 A 7^ 

10. / J 








1 7 9'^ 

i / .Z,J 






019 


1 7 7S 






01^ 
U.U i J 


1 S 9S 






O 01 ^ 


1 s 7^=; 

io. iJ 






01 s 
U.UI J 


19.25 


0.007 


0.025 


0.021 


19.75 


0.009 


0.006 


0.023 


20.25 


0.014 


0.010 


0.031 


20.75 


0.019 


0.013 


0.043 


21.25 


0.024 


0.021 


0.061 


21.75 


0.032 


0.028 


0.087 


22.25 


0.041 


0.031 


0.116 


22.75 


0.059 


0.041 


0.124 


23.25 


0.079 


0.056 


0.145 


23.75 


0.163 


0.078 


0.101 


24.25 


0.105 


0.087 




24.75 


0.111 


0.148 





tained by grouping the stars in magnitude bins and calculat- 
ing the differences Am = Wout ^ m-in between the retrieved 
and input magnitudes. The standard deviations are those of a 
Gaussian fitted to the observed distribution of Am. The Am 
are typically 50% larger than the DAOPHOT errors. The com- 
pleteness was calculated in a two-dimensional way by divid- 
ing each color-magnitude diagram in cells with size 0.2 mag 
in magnitude and 0.5 mag in color Contours of equal com- 
pleteness, smoothed over 3x3 cells, are superimposed to the 
color-magnitude diagrams. 



Observations of standard stars from Landolt (1992) were 
used to calibrate the photometry. A value of 0.030 s was as- 
sumed for the shutter delay time, which implies a correction 
of 0.003 mag for the standard stars. The raw magnitudes were 
normalized to 1 s exposure time and zero airmass: 



2.5 log(tcxp + At) -/caX 



(1) 



where map are the instrumental magnitudes measured in cir- 
cular apertures of radius R = 6'.'6 (close to the photoelectric 



aperture employed by Landolt 1992), A Hs the shutter delay 



and X is the airmass. The adopted mean extinction coefficients 
for La Silla are ks = 0.235, ky = 0.135, and kj ^ 0.048. 
A fit of the normalized instrumental magnit udes t o the magni- 
tudes of the standard stars in the Landolt' s ( 1992 ) fields gives 
the following relations: 

(2) 
(3) 
(4) 
(5) 



B = b' + 0.121 {B-V) + 24.777 

V = v' + 0.0425 {B-V) + 25.145 

V ^ v' + 0.0381 {V~I) + 25.145 
I = i' - 0.0253 {V-I) + 24.049, 



where the coefficients were derived from the standard stars 
observed on the second night. The r.m.s. scatter of the residuals 
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Fig. 2. Comparison of the present data with the photometry of 
van de Rydt et al. (1991). The differences in V (top panel), 
-B — (mid panel) and V—I (bottom panel) are plotted for stars 
brighter than V = 21. The dashed lines represent the median 
zero-point differences (this paper — VDK91) 

of the fit (0.007, 0.008, and 0.008 mag in B, V, and / respec- 
tively) was assumed to represent our calibration uncertainties. 

Before applying these calibrating relations to our photom- 
etry, the instrumental profile-fitting magnitudes (measured on 
the sum frames) were converted to the scale of the standard 
star measurements. To this purpose, a sample of bright isolated 
objects were selected in all individual images, all neighboring 
faint stars subtracted, and the cleaned frames were used to mea- 
sure magnitudes through circular digital apertures with radii 
'''ap = 6'.'6 (the same as for standard stars). Aperture correc- 
tions were then calculated for each image as the median of 5 
to 8 independent corrections in each filter. The r.m.s. scatter 
about the median zero point is 0.015 in B, 0.011 in V, and 
0.012 in /. This uncertainty includes not only the errors on 
aperture measurements and corrections, but also the effects of 
extinction variations during the nights. Following aperture cor- 
rection, calibrated magnitudes were obtained by applying the 
calibrating relations in the form of a linear system solved by the 
Kramer method (a procedure equivalent to the more usual iter- 
ative method). The total zero-point uncertainties, obtained by 
quadratic sum of the aperture correction and calibration errors, 
are 0.017, 0.014, and 0.014 mag in B, V, and / respectively. 

Figure ^ shows a comparison of our photometry with the 
results of van de Rydt et al. (1991) for stars brighter than 
V = 21. We obtain the following median magnitude and color 
differences (this study - VDK91): AV = -0.036 ± 0.083, 
A{B -V) = 0.091 ± 0.116, A{B - /) = 0.108 ± 0.151, and 
A{V — /) = 0.060 ± 0.101. The errors are standard deviations 



Fig. 3. The V, {B — V) color-magnitude diagram of the 
Phoenix dwarf galaxy, from our master photometric catalog. 
No selection has been applied to the data. Also shown are the 
10%, 30%, 50%, and 70% isocompleteness contours 

of the differences, following a 3cr rejection of the outliers. Van 
de Rydt et al. (1991 ) report the result of their comparison with 



the data of Ortolani & Gratton (1988). They found mean differ- 
ences, in the sense VDK91 -OG88,of AV^ = -0.014±0.080 
and A{B - V) ^ -0.082 ± 0.100. Thus our {B - V) colors 
are somewhat redder than those of VDK91 and consistent with 
the results of OG88. Our slightly brighter V magnitude scale 
is explained by the larger reference aperture (indeed, using the 
OG88 aperture we would have obtained a zero point dimmer 
by 0.03 ± 0.01 mag). 

3. The color-magnitude diagrams 

Figure |^ shows the V, {B ~ V) color-magnitude diagram of 
Phoenix. The most notable feature is our detection of the hori- 
zontal branch (HB) of Phoenix atV ^ 23.7, discernible up to 
{B — V) — 0.0. An extended blue tail, similar to that of metal- 
poor Galactic globular clusters, is not evident. The HB is found 
at a magnitude level where our CMD is significantly incomplete 
and photometric errors are quite severe, so deeper data are cer- 
tainly important. The morphology of the HB will be analyzed in 
more detail in Sect. ^. Figure |3] also displays a red giant branch 
(RGB) sharply cut at V ^ 20.5, and a blue sequence extend- 
ing to luminosities brighter than the tip of the RGB, represent- 
ing the young star population detected by OG88 and VDK91. 
This blue star sequence, partially overlapping the HB region, 
is also noticeable in other c-m diagrams. Figure Q shows the 
color-magnitude diagram of Phoenix using the wider {B — /) 
color baseline. The trend of the isocompleteness contours in 
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(B-l) (B-l) 

Fig. 4. The V, {B — I) color-magnitude diagram including all stars in the field of Phoenix {left panel). Superimposed are the 
10%, 30%, 50%, and 70% isocompleteness levels. A dotted box marks the region chosen to represent the blue sequence of young 
stars. The right panel shows a similar diagram for stars in the comparison field 



Fig. illustrates the relatively high incompleteness of /-band 
photometry in the case of blue stars. The CMD of Phoenix is 
compared with the similar CMD of the control field in the right 
panel of Fig. ^. A comparison of the two diagrams shows that 
the foreground and background counts do not represent a ma- 
jor source of contamination for the most interesting regions of 
the CMD. There are no foreground or background objects bluer 
than {B — I) = 0.5, so that we are confident that all stars bluer 
than this limit belong to a young stellar population in the dwarf 
galaxy. We have therefore defined a young star sample by se- 
lecting all stars in the CMD region having 19.5 < V < 23.5 
and B - / < 0.5. 

The spatial distribution of the blue star sample is plotted 
in Fig. H along with the distribution of stars in the red giant 
sample (defined below). The galaxy center, marked by a large 
cross, was defined using the mode of the marginal distribu- 
tions of star coordinates. The star formation sites in Phoenix 
are clearly delineated in this figure. The young stars are found 
in clumps in the central region of the galaxy, the most promi- 
nent clump being the well-known "association" described by 
Canterna & Flower ( |1977| ) and studied by OG88. The spatial 
distribution of the different stellar populations in Phoenix will 
be further discussed in Sect. ^ This information is used here to 
define a "red giant sample" by excluding all objects in the two 
most prominent star forming regions (approximated by the two 
rectangles in Fig. |]). 





• • 
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Fig. 5. Spatial distribution of the samples of red and blue stars 
discussed in the text (they are represented by small crosses 
and filled circles, respectively). Two small boxes are used to 
schematically delimit the star forming region. A large cross in- 
dicates the center of the galaxy. 
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Fig. 6. The V, {V — I) color-magnitude diagram of stars in 
the RGB sample of Phoenix, which excludes all stars within 
the star forming regions. Also shown are the fiducial red gi- 
ant branch sequence s of G alactic globular clusters from Da 
Costa & Armandroff (199C ), spanning a metallicity range from 
[Fe/H]= -2.2 to [Fe/H]= -0.7. A distance modulus (m - 
M)o — 23.1 was adopted for Phoenix 



Fig. 7. The color-magnitude diagram of Phoenix stars farther 
than 1'.3 from the galaxy center A horizontal branch extended 
to the blue is evident at y = 23.73 (this level is indicated by a 
horizontal line). Also shown are the fiducial sequences of RGB, 
HB, and AGB stars in the CMP o f the Galactic globular cluster 
M3 (from Buonanno et al. |1988[) 



Figure || presents the /, {V — I) color-magnitude diagram 
of the stars in the RGB sample, together with the fiducial red 
giant branch es of G alactic globular clusters from Da Costa & 
Armandroff (1990). From blue to red, the globular clusters are 
M15, NGC 6397, M2, NGC 6752, NGC 1851, and 47 Tuc, 
whose metallicities are [Fe/H] =-2.17, -1.91, -1.58, -1.54, - 
1 .29, and -0.7 1 dex, respectively. We note the extremely well 
defined RGB cutoff of Phoenix, suggesting the lack of a domi- 
nant intermediate age population. 

An expanded view of the horizontal branch of Phoenix is 
provided in Fig. ^ This diagram was obtained by picking only 
stars in the outer region of the galaxy (R > thus excluding 
all objects in the star formation regions. This selection proved 
to enhance the contrast between the HB and the RGB, partly 
because photometric error are smaller in the outer, less crowded 
regions, but possibly also because of an intrinsic difference in 
the radial distributions of the RGB and HB populations (we 
will return to this point in Sect. |]). Figure ^ clearly confirms 
the presence of a moderately blue horizontal branch in Phoenix 
at F 23.7. We note that the HB is hardly detected in the V, 
{V — /) diagram because most of the faint, blue HB stars fall 
below the detection threshold in the / band. The fiducial loci 
of the Galactic globular cluster M3 (which appears obviously 
more metal-rich than Phoenix: [Fe/H]= —1.57) are shown for 
comparison purposes, appropriately shifted to th e Ph oenix dis- 
tance as derived from the RGB cutoff (see Sect. 4.1 ). The data 



for M3 are from Buonanno et al. (1988), and the apparent dis- 
tance modulus and reddening from Harris (1997). 



4. Distance and metallicity 

The sample of red giant branch stars used to construct the 
Phoenix RGB luminosity function (LF) is a subset of the mas- 
ter red sample defined in the previous section. To reduce the 
contamination by young stars and field objects, we further se- 
lected a "2(7 RGB sample" which comprises all stars within 
±2(7 from the fiducial ridge line. The details of the method are 



given in Sect. 4.4 . The stars in the 2a red giant sample are used 



below to derive the fundamental parameters of Phoenix, dis- 
tance and metal abundance. 



4.1. Luminosity function and distance based on the RGB tip 

The distance to Phoenix was first estimated from the / absolute 
magnitude of the tip of the red giant branch (see Madore & 



Freedman 1995 for a discussion of the method and previous 
work). The cutoff was found at /tip = 19.09 ± 0.05. The error 
mainly reflects the uncertainty in locating the RGB tip, which 
is larger than the systematic error on our / magnitude scale. To 
assess the effects of crowding on the estimated RGB cutoff, we 
simulated the upper part of the RGB using artificial stars in the 
CMD region 1.3 < - / < 1.7, / > 19.0. After processing 
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Fig. 8. A simulation of the effects of crowding and photomet- 
ric errors for stars with magnitudes and colors typical of the 
observed HB. Artificial stars originally distributed inside the 
outlined box are spread over a larger color range after process- 
ing. The median of the input and output V magnitudes are rep- 
resented by the dashed and solid line, respectively 



the artificial data as the real RGB stars, the cutoff magnitude 
of the retrieved stars was 0.03 mag too bright, independently 
of the radial distance from the galaxy center. A correction for 
this small bias was included in our estimation of the distance 
modulus. 

The bolometric and / luminosities of the RGB tip, M™*^^ 
and MJ^*^^, and the / bolometric correction, were then esti- 
mated using the relations of Da Costa & Armandroff (1990). 
From the extinction-corrected color of red giants near the tip, 
{V - /)o,tip 1-48 ± 0.04, we obtain BC'i = 0.52 ± 0.01. 
The bolometric absolute magnitude of the tip was then in- 
ferred from our estimate of the mean metal abundance of 



Phoenix, [Fe/H] = -1.81 ± 0.10 (cf. Sect. The result 
was Af™*^^ — —3.47 ±0.02 (the error reflects the uncertainty 
on metal abundance), from which Af/^^^ = -3.99 ± 0.02 
was finally obtained. For the extinction and reddening correc- 
tions we used Ay = 3.1 Eb-v = 0.06 ± 0.06 and Ey-i = 
1.28 Eb-v = 0.03 ± 0. 03, a dopting Eb-v = 0.02 mag 
from Burstein & Heiles ( 1982 ), with a ±0.02 mag redden- 
ing uncertainty. The distance modulus to Phoenix determined 
from the extinction-corrected magnitude of the RGB tip is then 
(m — M)o = 23.04 ± 0.07, where the error includes the pho- 
tometric and reddening uncertainties, and the uncertainty on 
jyjTRGB jjjjg distance modulus is slightly shorter than that 
derived by VDK91 using the same method and reddening. 

4.2. Distance based on horizontal branch stars 

Our first detection of the HB of Phoenix allows to obtain an 
independent estimate of the distance to Phoenix. The mean 
magnitude of the horizontal branch was calculated as the me- 
dian V magnitude of the 553 HB stars in the CMD region 
0.0 < B - V < 0.6, 23.4 < V < 24.0. The median is 
Vhb — 23.73, with a r.m.s. scatter 0.17 mag, yielding a for- 
mal error on the mean of ^0.01 mag. The mean value T4jb 
was found to remain constant for different radial sub-samples. 



within the internal errors. The total uncertainty, including the 
systematic error of the V magnitude scale, is cry = 0.016 mag. 

Given the significant incompleteness of our photometry at 
the HB level, the effects of errors and incompleteness on the 
measured location of the HB need careful investigation. The 
results of our artificial star experiments were used to simulate 
the HB as a naiTow strip in the CMD (23.5 < V < 23.9, 
0.0 < B ~ V < 0.6). The distribution of the retrieved arti- 
ficial stars in the color-magnitude diagram (Fig. ^ was then 
analyzed exactly in the same way as the real HB data. The me- 
dian HB level of the retrieved stars was V = 23.65, i.e. the 
simulated HB appears slightly biased toward brighter magni- 
tudes due to the rapidly changing photometric completeness, 
though the effect is relatively modest. Taking into account this 
bias, we adopt Vhb = 23.78 ± 0.05 as the mean observed mag- 
nitude of the blue HB. 

Using this value for Vhb, we calculated the distance mod- 
ulus of Phoenix on the Lee et al. (1990) distance scale, using 
their relation for the absolute visual magnitude of RR Lyrae 
variables. 



M^^ 0.17 [Fe/H] +0.82 



(6) 



for a helium abundance of F = 0.23. We adopt this cal- 
ibration since it is al so the basis of th e RGB tip method (Da 
Costa & Armandroff |199C|; Lee et al. |1993|). On the Lee's et 



al. scale, which gives My^ — Aly^ = 0. 51 m ag for a metal- 



Hcity [Fe/H]= -1.81 ± 0.10 (cf. Sect. we estimate a 
distance modulus (to - M)o = 23.21 ± 0.08. Note that the 
distance error includes the photometric (statistical and zero- 
point) errors and the reddening uncertainty, but does not take 
into account the uncertainty of the adopted HB calibration and 
distance scale. For example, using the empirical calibration of 
the mean absolute magnitude of the HB in 8 M3 1 globular clus- 
ters by Fusi-Pecci et al. ( 1996 ), the distance modulus would be 
(to - M)q = 23.01 ± 0.10, i.e. brighter by about 0.2 mag. 

The distance based on the HB is only marginally larger than 
that determined from the tip of the RGB, and the agreement is 
good in view of the many sources of uncertainty. The average 
of the two determinations, (m — M)o = 23.12 ± 0.08, corre- 
sponding to 421 ± 16 kpc, is finally adopted as our best estimate 
of the distance to Phoenix. 

4.3. Mean abundance 



As in Paper I, the mean metal abundance of Phoenix was esti- 
mated by direct comparison of the red giant branch in the V , 
{V — I) color-magnitude diagram with the RGB ridge lines 
of the template globular clusters from Da Costa & Armandroff 
( |199C| ). We calculated the mean color differences 5{V — I)o 
between the data in the Phoenix red giant sample and the fidu- 
cial loci of the template clusters, using different choices for 
the luminosity range. The values of 5{V — I)o were corrected 
for extinction adopting Ey^i = 0.03 ± 0.03. We note that 
for magnitudes fainter than / ^ 20.5, the color distribution 
of the Phoenix RGB is skewed towards bluer colors, due to 
the presence of the unresolved AGB of the old population (a 
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Fig. 10. The {V — I) color dispersion of red giant stars in 
Phoenix in 0.5 mag bins (open squares), plotted together with 
the instrumental color scatter derived from simulated stars 
(filled triangles) 



Fig. 9. The mean color difference of the Phoenix red giant stars 
from the fiducial sequences of template globular clusters is 
plotted against the mean cluster metallicity. The error bars rep- 
resent the standard deviations of the residuals in {V — I). Also 
shown is a linear regression in the range —2.2 < [Fe/H] < 
— 1.3 (dotted line). The large cross identifies the mean abun- 
dance of Phoenix together with its error. 



is probably to be ascribed to the small shift in the calibrated 
{V - I) colors (cf. Sect. |). Using a wider magnitude range 
for calculation of the 6{V — /)o values, a slightly lower metal 
abundance would have been derived (e.g., [Fe/H]= —1.92 us- 
ing the upper 2 mag of the RGB). This reflects the bias in the 
giant branch color due to the presence of AGB stars. The dif- 
ference is comprised in the quoted uncertainties, though. 



similar trend was noticed in the CMD of And I by Da Costa 
et al. 1996). In Galactic globular clusters, AGB stars can be 
separated from red giant branch stars below an absolute mag- 
nitude Mv ~ —1.5. In Phoenix this precisely corresponds to 
V « 21.5 or / « 20.5. 

Figure ^ shows the mean color residuals S{V — /)o (calcu- 
lated in the range —4 < Mj < —3) against the metallicities of 
the Galactic clusters. A linear fit to the differences, excluding 
only the metal- rich cluster 47/Tuc, gives the relation 



[Fe/H] = -2.86 SiV - /)o - 1.81, 



(7) 



applicable in the range —2.2 < [Fe/H] < —1.3. This rela- 
tion implies that the interpolated metal abundance of Phoenix 
is [Fe/H] ~ —1.81 ± 0.10 dex. The abundance (la) error is 
derived from the total uncertainty on the mean {V — I) color 
using Eq. ^ The uncertainty on S{V — I)o includes the statisti- 
cal error (~ 0.01 mag), the systematic error due to calibration 
uncertainties (0.02 mag), and the reddening uncertainty. We do 
not include the uncertainties related to systematic differences 
between abundance scales for the Galactic globular clusters. 
Note that for metal-poor populations, however, there is reason- 
able agreement between different abundance scales (e.g., Car- 
re tta & Gratton|l997|). 



The mean abundance derived here is slightly higher than the 
value measured by van de Rydt et al. (1991). The difference 



4.4. Metallicity dispersion 



Table 3. Mean V — I and color dispersion along the RGB 



/ 


V -I 




av-i 








RGB 


simul. 


intrinsic 


19.25 


1.43 


0.090 


0.044 


0.079 


19.75 


1.33 


0.079 


0.049 


0.093 


20.25 


1.22 


0.079 


0.052 


0.059 


20.75 


1.11 


0.100 


0.052 


0.085 


21.25 


1.03 


0.150 


0.072 


0.137 


21.75 


0.95 


0.153 


0.087 


0.126 


22.25 


0.90 


0.183 


0.127 


0.132 


22.75 


0.80 


0.217 


0.163 


0.143 



Previous work has suggested the presence of an intrinsic 
color width in the red giant branch of Phoenix. We find confir- 
mation of this claim in our data. Table ^ gives our results for 
the mean color and measured RGB width at different luminosi- 
ties (columns 2 and 3). The mode and standard deviation of the 
{V — /) color distribution of red giants were measured in 0.5 
/-mag bins on a rectified version of the RGB, using a Gaus- 
sian fit with a 3(7 clipping to discard field objects and young 
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stars. The instrumental scatter derived from artificial star ex- 
periments is given in column 4. The color dispersion of simu- 
lated stars, distributed along the RGB fiducial locus, was mea- 
sured with the same robust approach used to derive the RGB 
widths. The observed color dispersions and instrumental er- 
rors are compared in Fig. Clearly, the RGB color scatter is 
significantly larger at any luminosity than expected from mea- 
surement errors alone. The intrinsic {V — I) dispersions cal- 
culated as the quadratic difference between the measured and 
instrumental scatter are given in col. 5 of Tab. ^. The average 
color dispersion in the 4 brightest magnitude bins (/ < 21.0) 
is 0.079 ± 0.008 (r.m.s. error on the mean). It is interesting to 
note that the values are roughly constant above 7 = 21. Below 
this magnitude, the presence of AGB stars evolving away from 
the horizontal branch precludes reliable analysis. 

This evidence for an intrinsic scatter has been indepen- 
dently confirmed by a direct compariso n of o ur photometry 
with the data of van de Rydt et al. (1991). In the range 
19.0 < I < 19.5, the color differences for the stars in com- 
mon yield a standard deviation of 0.072 mag. Under the sim- 
ple assumption that the two samples have comparable errors, 
we obtain an instrumental scatter of 0.05 1 mag, which implies 
0.074 mag for the intrinsic scatter, in good agreement with that 
obtained from crowding simulations. We thus conclude that an 
intrinsic color scatter is indeed present in Phoenix, of the order 
o'intriV — I) = 0.08±0.01 mag. If we are to explain this scatter 
as due to an abundance spread, it would correspond to a metal- 
licity dispersion ctifc/h] = 0-23 ± 0.03 dex (using Eq. Our 
data thus seem to indicate that the RGB stars in Phoenix span a 
(±lcr) range in metal abundance -2.04 < [Fe/H] < -1.58 
dex, i.e. a range of ^0.5 dex. This confirms the relatively 
modest abundance range found by OG88, while VDK91 give a 
somewhat larger metallicity spread. 

An abundance scatter similar to that found in the Phoenix 
dwarf is well established in many dwarf spheroidals, where 
it has been co nfirme d also by spectr oscopic observations 
(Suntzeff et al. |1993[ Da Costa |1998|). In dSph's, a range 



)3j Da Costa 

in metal abundance is indicative of an enrichment process in 
which multiple stellar generations have taken place from in- 
terstellar gas enriched by previous episodes. However, a range 
in age may also affect the RGB color dispersion. As an ex- 
ample, we have calculated the mean {V — I) color difference 
between two isochrones with ages of 15 and 5 Gyr respec- 
tively, using the models of Bertelli et al. (1994i. The shift 
was calculated in the luminosity range —4 < Mj < —3.5, 
i.e. the interval used for calculating cr[Fc/H] ^ for two different 
metallicities, Z = 0.0004 ([Fe/H]= -1.7) and Z = 0.004 
([Fe/H]= -0.7). We find a mean shift A{V - /) « 0.07 for 
the metal-poor isochrones and A{V — I) w 0.03 for the metal- 
rich ones (younger stars are bluer). The effect of an age mix 
on the {V — I) colors of the red giants is therefore not neg- 
ligible, particularly in the case of metal-poor isochrones, for 
which the color difference near the RGB tip mimics a metal- 
licity variation of ~ 0.2 dex. Thus, depending on the details of 
the star formation and enrichment history of each galaxy, the 
effects of a younger age and a higher metallicity may partly 



cancel out, in which case the abundance spread inferred from 
the color dispersion would underestimate the metallicity range. 
An extreme example of such a combined effect may be at work 
in the Carina dwarf, which has a very narrow RGB indicative 
of little chemical enrichment (of the order 0.2 dex) although it 
has had several generations of stars (Smecker-Hane et al. 1 994 ; 
Hurley-Keller et al. 1998). Understanding the observed color- 



magnitude diagrams of a complex stellar system requires real- 
istic modeling of both its star formation and chemical enrich- 
ment histories (best performed on high-resolution data), and 
independent information on metal abundances such as that pro- 
vided by spectroscopy of individual stars. 

5. The stellar content of Phoenix 

5.1. The old population: horizontal branch morphology 

We begin our analysis and discussion of the star content of 
Phoenix with a few comments on our HB detection. The 
moderately blue horizontal branch unquestionably implies that 
Phoenix harbors a sizable old stellar population. Due to the rel- 
atively large photometric errors, stars on the red HB cannot be 
disentangled from the RGB, yet we can obtain some informa- 
tion about the HB morphology looking at Fig. ^ The horizontal 
branch of Phoenix appears mostly populated on the red side and 
moderately extended to the blue (0.0 < B - V < 0.8), with 
a hint of the RR Lyrae gap at B ~ V ^ 0.4. Thus it appears 
remarkably similar to the HB types of the dwarf spheroidals 



Leo II (Demers & Irwin |1993[ Mighell & Rich [1996K And I 
(Da Costa et al. 119961), Draco (Grillmair et al' 119971), and 



1998t see also Da Costa |1998| ). Given 



Tucana (Seitzer et al. 
the mean metallicity of Phoenix, this HB morphology implies 
a mild second parameter effect, since old halo clusters with 
[Fe/H]'^ —1.8 have blue horizontal branches. Under the hy- 
pothesis that the HB morphology of Phoenix is mainly driven 
by age, the theoretical HB models of Lee et al. ( 1994 ), and in 
particular their HB-type versus metallicity diagram, would in- 
dicate for its old population an age ^2 — 3 Gyr younger than 
that of the old halo clusters. However, it appears most likely 
that the HB of Phoenix is due to a mixture of stars of different 
ages. 

Recent studies have shown that the blue HB population in 
some dSphSs (e.g.. And I and Fornax) is less centrally concen- 
trated than the red HB stars and the red giants, most likely due 
to a gradient in the mean age of the galaxy populations (e.g.. 
Da Costa et al. 



1996 



Stetson et al. 1998). To investigate the 
presence of a radial population gradient in Phoenix, we have 
compared the radial distribution of the blue horizontal branch 
stars (23.4 <V < 24.0, 0.0 < S - < 0.6) to that of red gi- 
ants just above the HB (22.8 <V < 23.4,0.5 < B-V < 1.1). 
The surface density profile of the HB stars appears to be sub- 
stantially more extended than that of red giants. A two-sided 
Kolmogorov-Smirnov test clearly indicates that the two spatial 
distributions differ at a 99% confidence level. Since there might 
be some concern that apparent gradients be induced by differ- 
ent crowding and completeness in the two magnitude and color 
ranges, we applied the same test to the radial distributions of 
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Fig. 11. AGB star candidates in Phoenix, selected as redder 
than {B — I) = 3.0 and brighter than / — 19.5 (squares). Open 
symbols identify possibly blended stars. The tw o carbon stars 
spectroscopically identified by Da Costa ( 1994 ) are indicated 
by their numbers in VDK9 1 . Crosses represent field contami- 
nation in the same region of the CMD 



simulated stars in the same CMD regions. The cumulative radial 
density profiles of the blue-HB and RGB artificial stars look 
very similar Formally, the null hypothesis (the two data sets 
are drawn from the same parent population) cannot be rejected 
at any significant confidence level. We therefore conclude that 
the detected gradient cannot be produced by instrumental ef- 
fects. 

This extended distribution of the blue-HB sample (i.e. 
of the oldest stars) can have different implications. A radial 
change in the horizontal branch morphology could be due ei- 
ther to a radial change in the mean age of the stellar popula- 
tions, or to a higher mean abundance in the central regions, 
or both. Clearly, it is not easy to decide between these alter- 
natives on the basis of our data (the radial dependence of the 
HB morphology will be better investigated using HST photom- 
etry). We just note that the observed central concentration of 



young and intermediate age stars (see Sect. 5.3) indicates that 



multiple star formation episodes occurred preferentially in the 
central regions. 

5.2. The intermediate age population: AGB stars 

Previous studies of Phoenix have found a few very red stars 
located above the red giant branch tip. These may be AGB stars 



belonging to an intermediate age component. Da Costa (1994) 



less luminous than the brightest carbon stars in the ~ 8 Gyr 
old SMC cluster Kron 3, and considerably fainter than those 
in the populous young (1-3 Gyr old) clusters in the LMC. This 
would indicate that Phoenix did not form many stars since 8-10 



Gyr ago (Da Costa 1998). This suggestion, however, is based 



on a very limited sample of C stars, and brighter intermediate 
age stars may be identified with a complete census of possible 
upper- AGB stars in this galaxy. 

The relatively large field investigated in this paper gives us 
the possibility to better assess the contribution of an intermedi- 
ate age population in Phoenix, and to obtain some constraints 
on the star formation history of this dwarf galaxy. To this pur- 
pose, we selected stars brighter and redder than the RGB tip 
(/ < 19.5, B — I> 3.0) as candidate upper- AGB stars younger 
of 10 Gyr (Fig. ^^). The lower magnitude limit, V ~ 19.5, 
was chosen slightly fainter of the RGB tip to account for pos- 
sible large magnitude errors due to blending. All selected stars 
were visually inspected, and those with elongated shapes or 
other hints of nearby companions were flagged as less reli- 
able, although further study (in particular spectroscopy) will 
be needed to confirm the nature of any candidate. Figure [TI| 
also shows the foreground and background objects found in the 
same color and magnitude range (crosses). The brightest stars 
are obviously in excess over field objects. 

Taking into account field contamination and incomplete- 
ness, we counted 1 1 AGB stars in the CMD region just above 
the RGB tip (18 < < 19), in the color range 3.0 < 
{B — I) < 4.4. This interval is equivalent to a range in bolo- 
metric luminosity between A/boi = —3.5 (the RGB tip, where 
AGB evolution terminates in metal-poor globular clusters) and 
A^boi = —4.5, this brighter end roughly corresponding to the 
maximum AGB luminosity reac hed b y 3 Gyr old stars (Fro- 
gel et al. 1990 ; Marigo et al. 1996 ). For comparison, 156 
RGB + AGB stars are counted within 1 mag below the RGB 
tip (19 < / < 20). 

The occurrence of upper- AGB stars in dSph's generally in- 
dicates the presence of a population significantly younger than 
that of galactic globular clusters, and this seems to be the case 
also for Phoenix. Alternative explanations for the stars brighter 
than the RGB tip seem unlikely. For instance, large ampli- 
tude long-period variable stars (LPV's) have not been found in 
old halo globular clusters as metal-poor as Phoenix (Frogel & 



Elias 1988; Frogel & Whitelock 1988). Even accounting for the 



modest abundance dispersion inferred in Sect. there is no 
evidence for a significant component in Phoenix with metal- 
licity higher than [Fe/H]« —1. We have also considered the 
possibility that some of our AGB candidat es are artifacts of 
blending of red giants near the tip (Renzini 1998). Assuming 
that the luminosity of Phoenix (A/y « —9.7, van de Rydt et al. 



obtained spectroscopic confirmation of two carbon stars with 
A/boi ~ —3.7 in Phoenix. These two confirmed C stars are 



1991 ) is uniformly distributed over a 4' x 4' area, and adopting 
a resolution element 1.7 □", we expect less than one resolution 
element (0.4) in our frame to contain two red giants near the 
RGB tip. Further, in our simulations we find just 1 star scat- 
tered above the RGB tip by photometric errors. 

Thus we conclude that most of the stars brighter than 
/ = 19 trace the extended AGB of an intermediate age pop- 
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ulation. Our counts of upper-AGB stars are used to estimate 
the contribution of the intermediate age component to the lu- 
minosity of Phoenix, following the methods of Renzini & Buz- 



zoni (1980 and Renzini (1998). The "fuel consumption the- 
orem" gives the number of stars in a post-MS evolutionary 
phase j, jij ~ B{t) Lt tj, as a function of the total bolo- 
metric luminosity of the population (Lt) and the lifetime of 
the phase (tj). The specific evolutionary fiux B(t) is of the 
order 2 x 10^^^ stars Lq'^ yr^^ in good approximation for 
any age between 3 a nd 10 Gyr. Even though evolutionary pre- 
dictions in Renzini ( 1998 ) are given for a stellar population 
with solar composition, the results for nj are almost unchanged 
for metal-poor systems because a lower metallicity has op- 
posite effects on the bolometric corrections and lifetimes of 
red giants at a given age (Maraston 1998, priv. comm.). For 
thermally pulsing AGE stars we assume a standard lifetime 
of 1 M yr for each magnitude of luminosity increase (Renzini 
1998 ). An absolute magnitude My ~ —9.7 is adopted for 
the red populations of Phoenix, corresponding to a total lu- 
minosity Lv = 6.2 X 10^ Lq (van de Rydt et al. |1991| ). If 
the bulk of the Phoenix population were of intermediate age 
(3 < t < 10 Gyr), then the expected ratio between the number 
of thermally pulsing AGB stars in the 1 mag interval above 
the tip, and the sum of RGB and early-AGB stars in the 1 
mag range below the tip, would be A^agb/^rgbt ~ 0.19. 
This is about 3 times the fraction observed in Phoenix, where 
^Vagb/A^rgbt = 11/156 = 0.07 ± 0.03 (the error simply re- 
flects the count statistics). The number of upper-AGB stars is 
therefore consistent with '--^ 37 ± 12% of the stellar popula- 
tion in Phoenix being of intermediate age, i.e. even the small 
observed number of AGB stars implies a significant intermedi- 
ate age population. This not surprising since the lifetime of the 
thermally pulsing AGB phase is short enough that a large sam- 
ple of progenitors is required to observe an appreciable number 
of upper-AGB stars. The possible concentration of stars on the 
red side of the HB is indeed consistent with the presence of an 
~ 8 Gyr component, although a firm quantitative estimate of 
the fraction of intermediate age stars, and their mean age, must 
await for deeper images. 

5.3. The young stars 

We now turn to one of the most distinctive features of Phoenix, 
the presence of a recent burst of star formation. The wide base- 
line of the {B — I) color used in this study, and the avail- 
ability of a comparison field, allow us a better separation of 
stars in different evolutionary phases. Figure |l2| presents our 
CMD for stars comprised within the star formation regions (the 
rectangles in Fig. |]). Superposed on the diagram are three rep- 
resentative isochrones from the Padua stellar evolution mod- 
els (Bertelli et al. |1994| ), for a metal abundance Z = 0.0004 
([Fe/H]« -1.7) and ages 1.0 x 10^ 2.5 x 10^ yr, and 6.3 x 10^ 
yr. The corresponding masses at the MS turnoff are 4.8, 3.0, 
and 2.0 Mq, respectively. The model colors are sensitive to the 
adopted metallicity particularly in their evolved part. A ~ 10* 
yr old model provides a good fit to the upper main sequence of 
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Fig. 12. The color-magnitude of Phoenix in the star-forming re- 
gions, compared to a set of theoretical isochrones from Bertelli 
et al. ( |1994| ), with metallicity Z = 0.0004 and ages (1.0, 2.5, 
6.3) xlO* yr The dashed lines represent the main sequence 
evolution. Also shown are the foreground/background objects 
in an equal area of the comparison field (crosses) 



the young population in Phoenix. While in general core helium 
burning (HeB) stars on the blue loops give a substantial con- 
tribution to the b lue plume in ground-based photometry (e.g., 
Tosi et al. |1991|), the gap at {B - I) « 0.5 in Fig. O sug- 



gests that we have reached a good separation of the hydrogen 
and HeB star sequences. The group of bright blue stars having 
[B — I) ^ 0.5 and V < 21.5 matches quite well the expected 
location of HeB stars belonging to a 250 Myr old population. 
This suggestion is confirmed by the presence of a distinct se- 
quence at [B — I) ^ 1.5 which appear to trace the red end 
of the blue loops. Also, a short sequence of stars comprised in 
the region 0.5 < [B - I) < 1.5, 22 < y < 23 is quite well 
matched by the blue loop phase of the 6 x 10* yr old isochrone. 
The brightest red stars iV < 20) almost overlap the red super- 
giant phase of the same 6 x 10* yr isochrone (although the 
precise location of these luminous stars is a sensitive function 
of metallicity). Thus, the emerging picture is one in which the 
most recent star formation episode in Phoenix had finite dura- 
tion (of the order 0.5 Gyr or more), with evidence of a distinct 
burst ^0.6 Gyr ago, and went on forming stars until about 10* 
yr ago. It is interesting to note that this picture is similar to what 
we have found in Fornax, where star formation proceeded from 
1 Gyr until 0.2 Gyr ago (Held et al. 



1999) 



How significant is this young population in terms of lumi- 
nosity and mass ? To estimate the V luminosity of the young 
star component, we selected the blue stars in the color and mag- 
nitude range -1< B - I < 0.5, 21.5 <V < 23.0 (the faint 
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The spatial distribution of young and intermediate age stars, 
and their relationship with the neutral gas in Phoenix (Young & 
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Fig. 13. Plot of the spatial distribution of the Phoenix young 
and intermediate age populations. Open circles represent blue 
stars, with larger symbols indicating stars brighter than V = 
22.5. Crosses are upper-AGB stars (a few big crosses indicate 
objects possibly blended). We also partially reproduce the (0.5, 
1, 2) xlO^^ cm~^ Hi column density contours from the maps 



ofYoung&Lo(1997) 



limit roughly corresponds to the MS turnoff of the 0.25 Gyr 
isochrone). The stars were counted in 0.5 mag bins, and cor- 
rected for incompleteness according to the results of the arti- 
ficial star experiments. We found 57 blue stars, corresponding 
to 106 stars after completeness correction. Assuming that all 
these are MS stars distributed according to a Salpeter initial 
mass function, an order-of-magnitude estimate of the V lumi- 
nosity of the youngest stars was obtained following Renzini 
(1998), under the simple assumption of a single 100 Myr old 
burst of star formation. The observed stars are predicted in the 
mass range 3.0 < M < 4.8 Mq for a total V luminosity of 
Lv ~ 3.6 X 10'' Lq. Using the bolometric corrections and 
mass-to-light ratios of Maraston (1998), these correspond to 
^ 5.6 X 10^ Lq in bolometric units or a mass of ^ 7 x IQ-^ 
Mq for the most recent burst. The youngest stars thus con- 
tribute <^6% to the V luminosity of the galaxy, corresponding 
to a mass fraction Myoung/Afoid ~ 0.002. These results are 
similar to those obtained by Mould (1997) and Aparicio et al. 



(1997) for LGS3, a dwarf galaxy which shares the basic char- 



acteristics of Phoenix, i.e. a smooth optical appearance accom- 
panied by recent star formation. These estimates do not include 
the blue stars older than 250 Myr, in particular the 600 million 
yr component. The star formation history of Phoenix (in partic- 
ular in the last 1 Gyr) will be studied in more detail in a future 
paper using deeper data and synthetic CMD's 



Lo 1997), are shown in Fig. |13[ The original blue star sample 
(Sect. ^ was split at — 22.5 in two subsets of stars older and 
younger than ~ 2 x 10^ yr (open circles; larger symbols repre- 
sent the brighter, younger stars). The location of the youngest 
stars suggests that the most recent star formation took place in 
a few central sites, while stars older than ^0.2 Gyr are spread 
over a larger area. The overall distribution of blue stars appears 
elongated in the NE-SW direction, i.e. roughly perpendicular 
to the halo of Phoenix (cf. VDK91 and Fig. |l]) and pointing to- 
wards the Hi cloud. Since the radial velocity of Phoenix is not 
known, we do not really know if the Hi cloud A is physically 
associated with the galaxy, or the Magellanic Stream, or neither 
(e.g.. Young & Lo 1997). However, if we assume that cloud A 
is associated with the galaxy (as its location and shell-like ap- 
pearance seem to suggest), then the mass in neutral gas would 
be 1.2 X 10^ Mq. This corresponds to Mni/iy « 0.04, an 
order of magnitude smaller than the Hi content in dwarf irreg- 
ulars, which typically have a gas-to-star mass ratio of the order 
of unity (Mjjj/Ly = 0.8 - 3; Carignan et al. |1998[ ). Among 



the several plausible scenarios, we only briefly comment on the 
possibility that the cloud A consists of gas accumulated from 
mass lost by evolved stars, blown out by the energy input of the 
most recent star formation episode. We have made a simple cal- 
culation of the mass released by an old population by adopting 
a total mass loss rate from red giants, AGB stars and plane- 
tary nebula ejecta of 0.015 Mq yr~^ per 10^ yr (cf. Mould et 
al. |1990| ). We obtain a gas return of ~ 7 x 10^ Mq after 10 
Gyr of normal evolution of old stars in Phoenix, i.e. the ob- 
served amount of Hi could have been accumulated in 2 Gyr. 
By comparison, we estimate that only ^ 10'^ Mq of gas have 
been returned to the interstellar medium by type II supernovae 
evolved from massive (> 8 Mq) stars in the recent burst. Al- 
though the ability of a dwarf galaxy to retain its gas depends 
on poorly known physical parameters, such as the effects of 
SN explosions and the presence of a dark matter halo, these 
figures are not inconsistent with an internal origin of the gas in 
Phoenix. 



6. Summary and conclusions 

We have presented a deep CCD study of the Phoenix dwarf, 
a galaxy often regarded as a transition case between gas-poor 
dwarf spheroidals and gas-rich dwarf irregulars. Here we sum- 
marize the main conclusions of this study. 

The detection of the HB of Phoenix represents one of the 
main results of this paper. We find a mean magnitude Vhb = 
23.78 ± 0.05 after correction for instrumental biases. Informa- 
tion on the HB morphology was obtained using a statistical ap- 
proach. The horizontal branch of Phoenix turns out to be well 
extended to the blue, although red stars are about twice as nu- 
merous. This morphology, similar to that of And I, Leo II, and 
Tucana has two important implications. First, it demonstrates 
the presence of a significant (if not dominant) population older 
than 10 Gyr. Second, it implies that Phoenix, given its low 
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metallicity, represents yet another (mild) example of the "sec- 
ond parameter effect" in dwarf galaxies. If the origin of this 
effect is identified with an age difference, the bulk of the stellar 
populations in Phoenix would have to be younger (by 2-3 Gyr) 
than stars in old halo globular clusters. We find that the spatial 
distribution of blue (B — V < 0.6) HB stars is significantly 
more extended than that of red giant stars. Our result indicates 
that the early star formation episode occurred in Phoenix on a 
larger spatial scale that subsequent bursts. Horizontal branch 
stars have been used as tracers of the spatial distribution of the 
oldest stellar populations in other dwarf galaxies (And I, For- 
nax, Antlia), with similar results. 

Besides this old population. Phoenix has a significant inter- 
mediate age component. We have confirmed the presence of a 
small number of stars above the RGB tip, significantly in ex- 
cess over field contamination, and argued that these most likely 
trace the extended AGB of an intermediate age population. Us- 
ing a standard lifetime for the upper AGB stars, we estimate 
that the intermediate age populations contribute about 30-40% 
of the V luminosity of Phoenix. While we cannot establish the 
age of this component, there is some evidence that star forma- 
tion declined since 8-10 Gyr ago. The candidate AGB stars 
seem to concentrate in the inner part of the galaxy, although to 
a lesser degree than the young stars. 

Our wide photometric baseline has provided new informa- 
tion on the young stellar population in Phoenix. We have shown 
that the recent star formation episode, responsible for the sprin- 
kling of blue stars characteristic of this dwarf, started at least 
0.6 Gyr ago. The recent burst of stars formation (~ 1—2.5x10^ 
yr ago) accounts for less than 6% of the V luminosity of 
Phoenix and 0.2% in terms of mass. The blue stars which trace 
the most recent burst are concentrated in clumps or "associa- 
tions" near the galaxy center, with a spatial distribution elon- 
gated in a direction perpendicular to the major axis defined by 
the diffuse galaxy light, and slightly offset towards the Hi cloud 
observed by Young & Lo (1997). The neutral gas could have 
been blown out by the recent burst, a possibility the should be 
further investigated when the hypothesis of a physical link re- 
ceives support by measurements of the galaxy radial velocity. 

Excluding the regions of recent star formation, we have de- 
fined a clean sample of RGB stars that has been employed to 
re-derive the galaxy basic properties. A new distance modulus 
(m — M)rgb,o = 23.04 ± 0.07 was obtained using the well 
defined cutoff of the red giant branch in the /, (V — I) diagram. 
More importantly, we have obtained for the first time an inde- 
pendent estimate of the distance to Phoenix, (m — Af )hb,o = 
23.21 ± 0.08, based on the mean level of horizontal branch 
stars. The mean of the two independent measurements gives a 
distance (m — A/)o = 23.12 ± 0.08, which confirms previous 
estimates. Amean metal abundance [Fe/H]= — 1.81±0.10 was 
obtained with a direct comparison of the upper part of the RGB 
with the fiducial sequences of template Galactic globular clus- 
ters. A careful analysis based on extensive artificial star tests 
and comparison with previous photometry confirms the pres- 
ence of an intrinsic color scatter in the red giants of Phoenix, 



'^{v-i) = 0-08 ± 0.01 mag, corresponding to a metallicity dis- 
persion cr[Fe/H] = 0.23 ± 0.03 dex. 

In conclusion, this paper provides new evidence that 
Phoenix has had an extended history of star formation. Its stel- 
lar populations appear fundamentally similar to those found in 
many dwarf spheroidal galaxies - even the young population 
has ages comparable to those of the youngest stars in Fornax 
(Beauchamp et al. |1995| , Stetson et al. |1998| , Saviane & Held 



1999; Held et al. 1999). Its Hi content (if any) is compara- 



ble to the amount of neutral gas found in Sculptor. In view 
of these results, we are inclined to regard Phoenix as a low- 
mass dwarf spheroidal seen in the middle of a star formation 
episode (during which gas is perhaps expelled). This similarity 
in the stellar populations of Phoenix and more luminous dSph's 
is noteworthy given the difference either in mass or in location 
with respect to the big Local Group spirals. There appears to 
be no obvious correlation between the timescale of star forma- 
tion and galaxy mass, although there is one between mass and 
production of heavy elements, as made evident by the well- 
know luminosity-metallicity correlation (e.g., Buonanno et al. 
1985 ). The effects of the environment are not clear either, al- 
though Phoenix seems to fit well the trend suggested by van den 
Bergh (1994) between the presence of young or intermediate 
age populations in dSph's and their distance from the Galaxy 
or M3 1 . The origin of the striking difference in the star forma- 
tion histories of Phoenix and Tucana (both are isolated Local 
Group dwarfs) remains rather puzzling, since these two dwarfs 
not only have comparable luminosities and metal abundances, 
but also a similar HB morphology, which probably indicates an 
old population formed nearly at the same epoch. 

Acknowledgements. It is a pleasure to thank C. Maraston and 
G. P. Bertelli for clarifying discussions on stellar population synthesis 
in metal-poor systems. We are grateful to C. Chiosi and S. Ortolani for 
careful reading of the original manuscript, and to the referee, Dr. De- 
mers, for useful comments and suggestions. G. Da Costa kindly pro- 
vided unpublished information on his spectroscopic observations of 
candidate AGB stars. I.S. acknowledges support of ANTARES, an as- 
trophysics network funded by the HCM programme of the European 
Community. Y. M. acknowledges support from the Italian Ministry of 
Foreign Affairs and the Dottorato di Ricerca program at the University 
of Padova. 



References 

Aparicio A., Gallart C, Bertelli G., 1997, AJ 1 14, 680 
Beauchamp D., Hardy E., Suntzeff N.B., Zinn R., 1995, AJ 109, 1628 
Bertelli G., Bressan A., Chiosi C, Fagotto R, Nasi E., 1994, A&AS, 
106, 275 

Buonanno R., Corsi C.E., Fusi Pecci R, Hardy E., Zinn R., 1985, A&A 
152, 65 

Buonanno R., Buzzoni A., Corsi C. E., Fusi Pecci R, Sandage A. 
R., 1988. In: Grindlay, J. E., Davis Philip A. G. (eds.) Proc. lAU 
Symp. 126, The Harlow- Shapley Symposium on Globular Clus- 
ter Systems in Galaxies. Kluwer, Dordrecht, p. 621 

Burstein D., Heiles C, 1982, AJ 87, 1165 

Cantema R., Flower P J., 1977, ApJ 212, L57 

Carignan C, Demers S., Cote S., 1991, ApJ 381, L13 



14 



Enrico V. Held et al.: Stellar populations in the Phoenix dwarf galaxy. 



Carignan C, Beaulieu S., Cote S., Demers S., Mateo M., 1998, AJ 
116, 1690 

CarrettaE., Gratton R. G., 1997, A&AS 121, 95 

Da Costa G.S., 1994. In: Meylan G., Prugniel R (eds) Proc. ESO 

Workshop 49, Dwarf Galaxies. ESO, Garching, p. 221 
Da Costa G.S., 1997. In: Amaboldi M., Da Costa G. S., Saha P (eds.) 

Proc. ASP Conf. Ser. 116, The Second Stromlo Symposium: The 

Nature of Elliptical Galaxies. ASP, San Francisco, p. 270 
Da Costa G.S., 1998. In: Aparicio A., Herrero A. (eds.) Proc. of the 

Vlllth Canary Islands Winter School, Stellar Astrophysics for the 

Local Group: A First Step to the Universe. Cambridge University 

Press, Cambridge, p. 351 
Da Costa G. S., Armandroff T. E., 1990, AJ 100, 162 
Da Costa G. S., Armandroff T. E., Caldwell N., Seitzer P, 1996, AJ 

112, 2576 

Demers S., Irwin M. J. 1993, MNRAS 261, 657 

Frogel J. A., Elias J. H., 1988, ApJ 324, 823 

Frogel J. A., Whitelock P, 1998, AJ 1 16, 754 

Frogel J. A, Mould J. R., Blanco V. M., 1990, ApJ 352, 96 

Fusi Pecci F, Buonanno R., Cacciari C, Corsi C. E., Djorgovski S. 

G., Federici L., Ferraro F. R., Parmeggiani G., Rich M. R., 1996, 

AJ 112, 1461 

Grillmair C. J., Mould J. R., Holtzman J. A., et al. (WFPC2 IDT), 

1997, AJ 115, 144 

Hai-ris W. E., 1998, Electronically Published 

Catalog of Galactic Globular Clusters, 

|ittp://physun.physics.mcmaster.ca/Globular.htm: 

Held E.V., Bertelli G., Saviane 1., 1999, in preparation 

Hurley-Keller D., Mateo M., Nemec J., 1998, AJ 115, 1840 

Landoh, A. U., 1992, AJ 104, 340. 

Lee M. G., Freedman W. L., Madore B. F, 1993, ApJ 417, 553 

Lee Y. W., Demarque P, Zinn R., 1990, ApJ 350, 155 

Lee Y. W., Demarque P, Zinn R., 1994, ApJ 423, 248 

Madore B. F, Freedman W. L., 1995, A J 109, 1645 

Maraston C, 1998, MNRAS 300, 872 

Maiigo P, Bressan A., Chiosi C, 1996, A&A 313, 545 

Mateo M., 1998, ARA&A 36, 435 

Mighell K. J., Rich R. M., 1996, AJ 1 1 1, 777 

Mould J. R., 1997, PASP 109, 125 

Mould J. R., Bothun G. D., Hall P J., Staveley-Smith L., Wright A. 

E., 1990, ApJ 362, L57 
Ortolani S., Gratton R.G, 1988, PASP 100, 1405 (OG88) 
Renzini A., 1998, AJ 115, 2459 

Renzini A., Buzzoni A., 1986. In: Chiosi C, Renzini A. (eds.) Spectral 

Evolution of Galaxies. Reidel, Dordrecht, p. 195 
Saviane 1., Held E.V., 1999, in preparation 
Saviane 1., Held E. V., Piotto G., 1996, A&A, 315, 40 (Paper I) 
Seitzer P, Lavery R. J., Suntzeff N. B., Walker A. R., Da Costa G. S., 

1998, in preparation 

Smecker-Hane T. A., Stetson P B, Hesser J. E., Lehnert M. D., 1994, 

AJ 108, 507 
Stetson P B., 1987, PASP 99, 191 

Stetson P B., Hesser J. E., Smecker-Hane T. A. 1998, PASP 1 10, 533 
Suntzeff N. B., Mateo M., Temdrup D. M., et al., 1993, ApJ 418, 208 
Tosi M., Greggio L., Mai'coni G., Focardi P, 1991, AJ 102, 951 
van den Bergh S., 1994, ApJ 428, 617 

van de Rydt F, Demers S., Kunkel W. W., 1991, AJ 102, 130 (VDK91) 
Young L. M., Lo K. Y, 1997, ApJ 490, 710 



